SPECTRAL THEORY AND NONLINEAR PDE: A SURVEY 
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This paper is not intended as a sweeping survey on the topic of spectral theory 
and PDE as this would exceed any reasonable bounds. Rather, we intend to outline 
some of the recent work by the author and Joachim Kricger on the topic of stable 
manifolds for unstable PDEs. More specifically, we intend to give some ideas of how 
to prove a recent result by Krieger and the author for the quintic wave equation 

Hip- ip b = in M 3 . 

This is the focusing, iJ 1 -critical nonlinear wave equation (NLW). There exist sta- 
tionary solutions 4>(x, a) of this equation given by: 

<f>{x,a) = (3a) '(1 + a\x\ 2 )~i, a > 

The question is what happens if the initial data of NLW is a small perturbation 
of 4>. To answer this question, we first linearize around <f>. The ansatz ip(x t t) = 
4>(x, a) + u(x, t) leads to the equation 

d tt u - Au + Vu = N(u, <p), 

with V = — 50 4 and a nonlinear term 

N(u, (j>) = 10uV + lOu 3 2 + 5u 4 + u 5 

Let H = — A + V. It is known that the spectrum of such a Schrodinger operator on 
the interval (0, oo) is purely absolutely continuous. Also, H has a simple negative 
eigenvalue, say —k 2 (the ground-state). Indeed, 

inf (Hf, J) < (H(f), 0) = -4 / 6 dx < 

The presence of negative spectrum implies that the static solutions (/>(■, a) are lin- 
early (exponentially) unstable. We remark that by Agmon's estimate, we know 
that the ground-state function g, which satisfies Hg = —k 2 g, decays exponentially. 
In addition, g > 0, and g is radial. Because of symmetries, H has a zero mode. 
More precisely, by the dilation and translation symmetries 

Hd a <f>(;a) = 0, HV(p(-,a) = 

It is not an easy task to show that H has exactly one negative eigenvalue (which 
then must be simple), and no other zero modes than those four which we have 
specified. Numerically, this has been verified by means of the Birman-Schwingcr 
method (the numerics show that the Birman-Schwinger operator has exactly five 
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eigenvalues > 1 counted with multiplicity. This corresponds precisely to the four 
independent eigenfunctions and the resonance function). 

However, if we restrict H to the invariant subspace of radial functions, then the 
determination of the spectrum of this restricted operator is elementary. Indeed, 
Hf = Ef for radial / = f(r) is the same as (-£7 + V)rf(r) = Erf(r). In 
other words, we are reduced to an eigenvalue problem for the half-line operator 
H = — + V with a Dirichlet condition at r = 0. In our case, the function 
rd r 4>{r,a) has a unique positive zero and solves H(rd r (j)(r,a)) = (moreover, it is 
bounded and therefore a resonance function for H). By Sturm theory, it follows 
that there is a unique negative eigenvalue which is simple. It agrees with —k 2 
from above. Hence, the operator H has a one-dimensional exponential instability 
(in the radial setting) and the evolution of the wave equation driven by H has a 
codimcnsion one stable manifold. Our main theorem states that this remains true 
(in some sense) on the nonlinear level. 

THEOREM 1 

There exists a 5 > such that the following statement holds. Let Bs(0) be the 
ball centered at of size S in H^ ad (M. 3 ) x H 2 ad (M. 3 ). There exists a codimension 1 
Lipschitz manifold E which divides Bs(0) \ X into two connected components and 
such that if (u\, U2) G S n Bs(0), then the Cauchy problem for NLW with data 

(V(-,o),^(-,o)) = (0(-,i) + Ul ,« 2 ) 

has a global solution tp(x,t). Moreover, there is the representation 

tp(x, t) — (j>(x, a(oo)) + u{x, t) 
where |a(oo) — 1| < 5, and ||u('g£)||oo ^ ^(t)^ 1 as well as 

(u, d t u) = (v, d t v) + o i j lxL2 (l) as f^oo 

with Ov = and (v,d t v)(0) £ H 1 x L 2 . 

One motivation for this problem came from a numerics paper by Bizon, Chmaj, 
and Tabor which appeared in Nonlinearity in 2004. These are physicists who are 
mainly interested in the Einstein equations of general relativity, but proposed the 
quintic NLW as a model case. They predicted numerically that the stable manifold 
should exist, at least in a small neighborhood of (j>. Moreover, Bizon et. al. made 
conjectures about the behavior above and below the manifold (one should lead to 
blow-up, the other to dispersion). 

Another motivation for studying the quintic NLW came from the author's work 
on the focusing nonlinear Schrodingcr equation (NLS), which we now describe. 
Consider the equation 

id t V + A* + |*| 2ff * = in K 3 , 

for a < 2 (where 2 is the H 1 critical exponent in K 3 . We also have a C rit = | = §, 
which is the L 2 critical exponent). The Hamiltonian and the charge are conserved 
functionals: 
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W(*(t)) = 

Q(*W) = 




For cr < §, we get global existence in H 1 using the Gagliardo-Nirenberg inequality 
and the above conservation laws. For a > -j, we in fact can have blow-up in finite 
time. This result is easy to see simply by looking at the so-called virial identity. 
This approach is known as Glassey's argument and is based on the functional 

Indeed, if a > 2/d, then 

^ J \x\ 2 mt,x)\*dx<im(*(p)) 

Thus, any solution of negative energy must blow up in finite time (at least in the 
sense of the H 1 norm). Indeed, at some point the positive quantity above disappears 
and the solution must blow-up. Strictly speaking, to apply the virial identity, we 
need to work in the space (x) -1 !! 1 (it suffices to assume that *(0) € (x)^ 1 H 1 ). 

There exist stationary state solutions for NLS as well. Make the ansatz fy(t, x) = 
e ita ^ Then </> satisfies the equation 

(a 2 - A)0 = |0| 2 > 

Berestycki and Lions showed using variational methods that there exist solutions to 
such equations that are positive, radial, and exponentially decaying. In many cases, 
uniqueness of said ground-state solutions has been proved by Coffman, Kwong, 
and McLeod. </> must satisfy the Derrick-Pohozaev identities and for monomial 
nonlinearities we have a scaling property, namely 4>{x, a) = a~ 4>(ax, 1). 

There is a group of natural symmetries acting on NLS. There is the modulation 
symmetry \I/ — > e l7 \£' and the Galilei transforms 

Q(t) = e i (~ v2t + x - v ) e ~ i P-(v+ 2tv ) j 
for p = —i-§-- These are the unique transformations that satisfy 

G(t)e ltA - e itA £(0) 

with G(0) = e lx ' v e~' lp ' y . Note that this latter transform at time t — is the Galilei 
transform from classical mechanics, i.e., x h- > x + y, p i— » p + v. Finally, we have 
dilation symmetries from the scaling relation above for ground state solutions. In 
total, we have a 2d + 2 vector of symmetry parameters 

Ti" = (l,v,y,a). 

The relevance of these symmetries can be seen in the following theorem of Michael 
Wcinstcin. 



THEOREM 2 

For any e, there exists a 8 such that if 



inf||*(.,0)-e^(-+ 2 /,a)|| Hl < 8, 
i-v 
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then 



inf ||*(x,t) 



- e i7 0(- +y,a)\\ H i < e. 



It is essential to mod out the symmetries in the second inequality, since a small 
perturbation of the ground-state <fi (which unperturbed has zero momentum) will 
impart a nonzero (but very small) momentum on it. Consequently the soliton will 
have moved a large distance in finite (but large) time. In a similar fashion, the 
phase will also change by size one in finite time. 

Weinstein's theorem expresses what is known as orbital stability of the ground- 
state solutions. During the late 1980s Grillakis, Shatah, and Strauss developed a 
fundamental theory of orbital stability of standing waves generated by symmetries. 
Their theory applies to general Hamiltonian equations which are invariant under 
the action of Lie groups. In particular, it gives the dichotomy a < 2/d and a > 2/d 
for the focusing NLS, but applies to a much wider range of examples. 

Note that the conjectured instability of excited state solutions is largely open. 
An excited state is a solution of 



which changes sign. Important work has been done in that direction by T. P. Tsai 
and H. T. Yau for NLS with a potential. Generally speaking, NLS with a potential 
or Hartree equations is a fascinating field to which we cannot do justice here, see 
for example the work of Frohlich ct. al. 

Bcrcstycki and Cazenave showed that for a > § , blow-up can occur for arbitrar- 
ily small perturbations of ground states. This is easy for a = § using Glassey's 
argument from above and the fact that 7i((l + e)<f>) < for e > since Ti{4>) = 0. 

Now, let us discuss the spectral analysis of linearized operators for NLS. Set 
ip = e 4> + R and linearize. Note that we are assuming a does not change with 
time here, but we do not let that concern us at this point. We end up with an 
equation 



If 4> = 0, we have a self-adjoint operator with purely absolutely continuous spectrum 
given by (— oo, — a 2 } U [a 2 ,oo). If we do turn on 0, then this set turns into the 
essential spectrum of H by Weyl's criterion. In addition, there is an eigenvalue at 
due to the symmetries. Indeed, writing NLS as N^^) = where indicates the 
symmetries given by the vector tt acting on 'J, we conclude that TV' n = 0. 
Since clearly H^^) = Af 1 ^^) by definition of the linearized operator, it follows 
that d^Tr is a zero energy state of H. In the NLS case, they are (generalized) 
eigenfunctions with zero eigenvalue (but for the quintic NLW from above, there is 
a resonance and not just an eigenvalue at zero). If we look at R = u + iv and write 
H in terms of u, v, then we obtain the equivalent matrix operator 



A0 + a 2 0= Ic/f^ 




for N a nonlinear term quadratic in R and 
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where 

L_ = -A + a 2 - 2ff 

and 

L+ = -A + a 2 - (2<T + 1)0 2<T . 
Clearly, L-((fi) = and L + (V(/>) = 0. Hence, we have 

\. f V0 




G Ker (H) 



Differentiating L-<fi = in a yields L+(d a (/>) — —2a<f>. Therefore, 

d a <t> \ = ( 
I \ 2im 



H 



In fact, it is easy to see that 

\ f d a 




G Ker (H 2 



This shows that the eigenvalue is of geometric multiplicity at least d + 1 and of 
algebraic multiplicity at least 2d + 2, but we actually have equality here (at least 
if a 7^ 1) due to a result of M. Weinstein from 1986. So, Zf is truly a non-self 
adjoint operator with a nontrivial Jordan form at zero energy. In the critical case 
cr=| the dimension of the root space at zero increases to 2d + 4 due to the extra 
pseudo-conformal symmetry. 

Let us now analyze the operator e ltH . Using the L_, L + version of H, we see 
that Hf = Ef implies 

H 2 ( * W L " L + ° \(h\ =E *ffi 



h J \ L+L_ y V /a / \f* 
Since L-(f> = and (f> > it follows that L_ is a non-negative operator with as 
a ground state. By work of M. Weinstein, L + has exactly one negative eigenvalue 
and the kernel is spanned by the derivatives dj<f>. From the above argument, since 
E 2 fi is in the range of Z/_, we know that /i _L <fi. Therefore, we have 

y/LlL+y/Llf! = E 2 h, 



for fx = LJ 2 fx- Thus, E 2 is real as an eigenvalue of a self-adjoint operator. The 
dichotomy of linear stability or instability is decided by the question: When is 
E 2 < 0? In order to answer this, look at 



lnf (IT) < 

& inf (L+fJ) < 0. 
/j_0,||/|| 2 =i 

Denote this latter infimum as /j,q and assume that /iq < 0. One can show that the 
minimum is attained at some function /o. Using Lagrange multipliers, we have 

L +fa = Mo/o + A0, /o -L <£, Mo = (L+fo,fo)- 

If A = 0, then /io is a negative eigenvalue, hence the ground state of L + and this 
contradicts the orthogonality relation. So, we have fo = {L + — /.io) _1 0. Consider 
the function 

h(p) = ((L + -[i)- 1 q>^). 
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Note that this goes to — oo for \x tending to the unique negative eigenvalue of 
and h(0) makes sense since the kernel of L+ is orthogonal to 4> (remember it is V0). 

2 
2- 



We know h'(y) > as it is ||(L+ - /i)" 1 ^!!!- Hence, h(0) > iff /z < 0. But, 



h(0) = -^{8^ 



2a 
1 

4a' 

2 



Finally, one checks from scaling that <9 Q ||\I/ Q || < <s=>- a > J. This (up to 
the equality sign) is the well-known instability condition for NLS. We now state 
our first stable manifold theorem for focusing, cubic NLS in three dimensions. It 
is conditional on the assumption that there are no imbedded eigenvalues in the 
essential spectrum of the linearized operator. 



THEOREM 3 

(Schlag - 2004) 

Consider NLS with a = 1 in M 3 . Assume that the linearized operator does not 
have any imbedded eigenvalues in its essential spectrum. Then there exists S > 
such that the following holds: Let Bs(0) C W ' n W 1,2 be the ball centered at 
of radius S. There exists a codimension one Lipschitz graph £ C Bs(0) which 
divides Bg(0) \ £ into two connected components and so that for initial data in 
^(-, 0) G £ + (/>(■, 1), there exists a global solution 

tf(as,i) = W(x,t) + R(x,t), 

with ||i?(-,t)|| OCl <S(t)-i, R = e ltA R + o L 2(l). We have 

W(x,t)=e ie ^U(x-y(t),a(t)), 

6{t,x) = - ( v 2 {s)dx + v(t)x +7(t) 
Jo 

and 

y(t) = 2 f v(s)ds + y (t). 
Jo 

Finally, the vector ir(t) := (j(t) , v(t) , yo(t) , ct(t)) approaches a terminal vector 7r(oo) 
which is S close to the initial vector 7r(0) = (0, 0, 0, 1). 

It is clear that we cannot have stability for initial data in the entire ball <p+Bs(0). 
In fact, a = 1 is in the supercritical range for which Bcrestycki and Cazcnave have 
exhibited arbitrarily small perturbations of 4> that lead to finite-time blow-up. On 
the other hand, we do not understand the general behavior of solutions with data 
in^ + S tf (0)\E. 

We remark that the method of proof does not extend to the entire supercritical 
range a > 2/3. More precisely, an important gap property for the operators L± 
(which we defined above) breaks down below some value o% — 0.91396 . . .. The 
gap property here refers to L + and L_ having neither an eigenvalue in (0, a 2 } nor 
a resonance at the edge a 2 . This was shown in some recent numerical work by 
Laurent Demanet and the author. This numerical work also establishes that the 
gap property does hold above a* , and thus in particular for a = 1 (which is essential 
for the proof of the theorem). 
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For the proof of the theorem, we follow the basic strategy of the modulation the- 
ory of Soffer and Weinstein. More precisely, we seek an ODE for w(t) = (7, y, u, a)(i) 
and a PDE for R. The ODE is constructed in such a way that the non-dispersive 
nature of the root-space of H is eliminated. Note that the dimension of the root- 
space is 8, as explained above, and w(t) has eight components. This agreement is 
of course not only essential but also no coincidence, as 7r has as many components 
as there are symmetries, and the dimension of the root-space must also agree with 
the number of these symmetries. 

Of fundamental importance are the following estimates on the linearized opera- 
tor: 

||e itH P s /||oc < H||/||i, 

\\e ltH P s fh < H/lla 

where P s projects onto the stable part of the spectrum (i.e., it is the identity minus 
the Riesz projection onto the discrete spectrum). The proof of these bounds not 
only requires that there arc no imbedded eigenvalues in the essential spectrum, but 
also the following properties: 

• zero is the only point in the discrete spectrum 

• the edges ±a 2 are neither eigenvalues nor resonances 

Wc prove these facts by modifying some arguments of Galina Perelman for the 
case d = 1 and a = 2. It is at this point that the aforementioned numerical work 
by Demanet and the author becomes essential. 

In one dimension a corresponding result was proven by Joachim Krieger and 
the author in the entire supercritical range a > 2. This result was completely 
unconditional and also did not rely on any numerics. Rather, one can prove all 
the necessary spectral properties by means of the scattering theory established by 
Buslaev and Perelman, as well as by invoking the explicit form of the ground state 
in d = 1 (the latter is needed to prove that there are no imbedded eigenvalues). 
The linear estimates needed for d = 1 are 

\\e itH P s f Woo <t-?\\fh 
IK^-V^P./IU <t-t||(a:)/||i 

The latter bound is somewhat curious, in so far as it does not hold in the free 
case H = —d 2 . Rather, it exploits that the edges of the essential spectrum are not 
resonances (which of course fails in the free case). Heuristically, one can arrive at 
the second bound as follows. Suppose the scalar operator —d 2 + V does not have a 
resonance at zero energy. This is well-known to be equivalent to the fact that the 
distorted Fourier basis {e(x, C)}?gk vanishes when £ = 0. Now write the evolution 
relative to this basis: 

<'"' /•-/•: <•:. J e^ 2 e(x,0(f,<;0)dC- 

By the vanishing of e(-,0), we integrate by parts in £ gaining a factor of t . 
However, this evidently comes at the expense of a factor of x. Hence, we obtain the 
decay by invoking the dispersive decay on top of the i _1 gain, leading to the 
above bound. On a more heuristic level, we remark that compact data / should be 
thought of as having traveled distance t after time t. Due to the weight (x)^ 1 this 
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will lead to a gain of t^ 1 on top of the dispersive decay of The point of not 
having a resonance at the edge is to avoid having some mass " being stuck" around 
the origin. This would of course invalidate the heuristic idea of transport. 

The possibility of obtaining improved decay bounds at the expense of weights was 
apparently discovered by Murata in the 1980s. It applies only to those cases where 
the free operator has a zero energy resonance, i.e., when d = 1,2. In dimension 
two, Murata discovered an improvement of (logi) -2 over the free dispersive decay. 

We remark that the presence of these weights leads to substantial technical 
problems in the proof of the 1-d stable manifold theorem (as compared to the 3-d 
result). 

Finally, we would like to point the reader to the second stable manifold paper 
with Krieger which deals with the pscudo conformal blow-up solutions for critical 

3 

NLS in d = 1. The bound again plays a crucial role, but this is not the venue 
to explain the blow-up argument in any detail. 

As stated above, we desire to understand stability of the focusing, quintic NLW 

Uip - ip b = 

around the special solutions <p(-,a). Recall that Aubin showed that there exists 
4> = a) > such that —A<p(; a) — 5 (-, a) = 0. Explicitly, with a > 

<f>{x,a) = (3a)^(l + a\x\ 2 )-^ 
Note from the formula above that <fi(x,a) ~ t^t and also d a (f>(x,a) ^ jjj- In 

particular, neither <\> nor d a 4> belong to L 2 (R. 3 ). We also remark that it is essential 
for us to work with the quintic, i.e., if ^critical equation for a number of reasons. 
For example, by a result of Gidas and Spruck, there are no positive solutions of 
A0 + (f)P = in R 3 if p < 5. 

We wish to study radial perturbations to the initial data of the form 

^(0) \ _ f <t>(; l ) \ , f h 

M(o) )-{o J + \f2 

As mentioned before, Bizon, Chmaj, and Tabor, as well as Szpak, studied this 
problem numerically. They tested initial data (/i, /ij) from a one-parameter family 
of functions (e.g., Gaussians with changing variance) and observed that there is a 
unique point on such a curve of data at which there is long-time stability. This 
point splits the curve into two connected components (i.e., intervals). One interval 
led to instability (or possibly blow-up), whereas another led to scattering to a free 
wave. By letting this curve of initial data vary, they then observed that the unique 
stability point sweeps out a codimension one stable manifold on which the family 
4>{-,a) acts as a one-dimensional attractor. 

Our theorem establishes the existence of such a codimension one stable manifold. 
However, at this point we cannot describe the behavior off the manifold. We now 
state the theorem, which already appeared above, in full detail. 



THEOREM 4 

(Krieger-Schlag - 2005) 
Let R > 1. Denote 



X R = {(/i,/ a ) G Hl ad x H 2 radl suppifj) c B(0,R)}. 
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Define 

So = {(/i, /a) G : (*(•> l)/i + / 2 ,.9(-, = 0}, 
where ~k 2 is the unique negative eigenvalue of the linearized Hamiltonian and g is 
the associated ground-state eigenf unction. There exists S > such that if (/i, / 2 ) G 
■65(0) C So, then there is a real number /i(/i,/ 2 ) such that 

\Kiuh)\ < W(h,hWx R , 

\Kh,h)-h{h,h)\ < 5||(/i,/a)-(/i,/2)||j 



/or(/i,/ 2 ),(A,/ 2 ) e S 5 (0) 



□V> - -0 s = 

0(-,o) \ _ / 1) \ , / /i(-) + M/i,/2M-,i) 



and 



atvo.o) y v / V /at) 

ftas a global solution. Moreover, this global solution admits the representation 
ip(-, t) = (/)(-, a(oo)) +«(•, i) where 

H;t)\\oo < 6(t)-\ 
\\(v,d t v)-(v,dtv)\\ 6 i xLa -> 



t — > 00 where \3v = and (u, 9^5) g i/ 1 x L 2 . 



The radial assumption is more than a mere convenience. It assures that the 
bulk term remains at the origin. In other words, in case of nonradial perturba- 
tions we would need to allow (f> to move away from the origin, i.e., (j>(X(t, x), a{t)) 
where X(x,t) preserves the Lorentz invariance. In the same spirit, recall that the 
cigcnspacc at zero is generated by the functions V0 (an expression of translation 
invariance), whereas the resonance at zero results from the dilation invariance. The 
latter is the only allowed symmetry in the radial case. 

The map 

fl,(0) - # r 3 ad x ff r 2 ad : (A, / 2 ) 1 ^ + h(h, f 2 )g(; 1), f 2 ) 

traces out a (small) manifold E. By our estimates, E can be written as a Lips- 
chitz graph and Eo is the tangent plane to the desired manifold E (because of the 
quadratic bound on h). 

For a general Schrodinger operator H = — A + V in R 3 the spectrum changes 
under the scaling V — > X 2 V(Xx) by multiplication by A 2 . By inspection, <j>(\x, a) = 
\~^4>{x, \ 2 a). Hence, our linearized operators 

H{a) = -A-50 4 (-,a) 

have the property that the spectrum depends on a via scaling. In particular, this 
is true of k (and of course (?(•, a) depends on a in a similar fashion). 

To prove the theorem, we make the ansatz ijj(t,x) = <p{x,a(t)) + u(t,x) with 
some path that satisfies a(0) = 1 (we will in fact also need that d(0) = as well as 
\a(t)\ < 5(t)~ 2 ). Then we plug this ansatz into our NLW which yields 

dttu -Au- 5(/) 4 {-,a{t))u(t) = -d U (f>{-,a{t)) + N ((/)(■, a(t)) , u) 

with the nonlincarity 

N(u, <t>) = lOu 2 3 + lOu 3 2 + 5u 4 + u 5 
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Let us define V(a) = — 5</> 4 (-, a). We want to eliminate the time dependence in V, 
so we look at V(a(oo)) and define 

□v(a(oo))" = d tt u - Au - 5<f)(-, a(oo))u. 

Thus, the linearized equation becomes 

nv( a(oo ))U = -0tt0(-, o(t)) + o(t)), u) + {V(; 0(00)) - !/(-, a(t)))u =: (*) 

Before we say more on how to analyze this equation, let us look at the Schrodingcr 
operator H(a) = —A + V(-, a). The following relations were already noted above: 

H{a)d a (j) = from dilation invariancc 
H(a)V<f) = from translation invariance. 

Since d a <t> & L 2 we have a resonance at the point of the spectrum. 

Let us discuss the spectrum of the above operator (on the radial subspace). 
First of all, since V(a) decays like |s| -4 , a result of Agmon and Kato implies that 
there is no singular continuous spectrum as well as no imbedded eigenvalues in the 
continuous spectrum. We also know that the absolutely continuous spectrum is the 
same as that of —A, namely [0, 00). Hence, we can decompose L 2 . ad = L 2 p © L 2 C . 
The pure-point part L 2 p is a one-dimensional space consisting of spanjg}. Since g is 
a ground-state solution, it decays exponentially. This decomposition is blind to the 
resonance. In particular, there is no L 2 projection associated with the resonance 
function. We shall say more about the characterization of resonances in terms of 
the perturbed resolvent in the next lecture. The basic principle is as follows: Look 
at (H(a) — z 2 )" 1 . As z — > in Imz > 0, this can either approach a bounded or an 
unbounded operator in a suitable weighted L 2 topology. Due to the work of Jensen 
and Kato from 1978, we know that this operator can be Laurent expanded in the 
form 

{H (a)-z 2 )- 1 = 2 ^ + — + 0(1). 
z z 

Clearly, — ci must equal the orthogonal projection onto the zero energy eigenspace, 
whereas c_i has contributions from both the resonance and the eigenfunctions at 
zero. For now, let us not concern ourselves with respect to which topology this is 
valid. Let us rather note that in our application we not only eliminate c_2 but 
also that c_i is solely due to the zero energy resonance (recall the restriction to the 
radial subspace). 

Wc now return to a more detailed description of the main ansatz. Rewrite the 
linearized equation as a Hamiltonian system: 

d t U = JH 00 U + 

where 




H(a(oo)) 
1 
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It is of course necessary to understand the propagator e tjn °° , and therefore also 
the spectrum of the operator JTioo- To this end, note that if 

then ff(a(oo))A = -£ 2 /i- 



The conclusion here is that spec(J7ioo) = y/—spccH(a(oo)) = iRU {±fc}- More- 
over, both ±fc are simple eigenvalues with eigenfunctions G±. 

Define P_ and P + to be the Riesz projections associated with the eigenvalues 
— k and k, respectively. This can be done in the standard way by taking contour 
integrals of the resolvent. More importantly, they turn out to be very explicit rank 
one operators with kernels defined by the ground state function and eigenvalue. 
Decompose the full solution U(-, t) as 

U(-,t) = n+(t)G+ + n_(t)G_ + U(-,i), 

where (7 — P+ — P-)U = U. Then, we get the following hyperbolic system of ODEs 
in the plane (with k = fc(-, a(oo))): 

h + (t) \ ( ~k \ / n + (t) \ _ / F + (t) \ 

Now let us find the unique line in the plane with the property that initial data 
chosen on it lead to bounded solutions: 



n+(t) = e tk n + (0)+ / e^-^ F+(s) dx 
Jo 

(1) e- tk n + (t) -> = n+(0)+ / e~ ks F + (s)ds 

Jq 

/>OC 

=>n+(t) = - e k{t - s) F + (s)ds. 



It follows from the final relation that n + (t) inherits the decay of F + (t) (e.g. any 
decay like t~P with (3 > 0) provided n+(0) is chosen as in QJ. Thus, is precisely 
the stability condition that we use to define the (unique) correction to the initial 
data which leads to bounded U at each step of the Banach iteration. Summing up 
all these corrections then yields the function h(f\, A). Next, by the explicit form 
of the projections P± one checks that 



U 



u 

d t u 



where u = F , , „m. Hence, we can go back from the Hamiltonian system to the 
wave equation 

V (.,a(oo))U = ^. Q(oo)) (*), 

«(0) = P g i,a(oo))(fl + Kh,h)9(;l)), 

dtu(0) = P 9 |, o(oo)) /2 

where (*) stands for the right-hand side of the linearized equation, sec above. As 
a starting point for bounding the solutions of this wave-equation, we show that for 
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some constant Co ^ 

(2) sm (V g ( fl M)) f i / = Co5q 0(., a(oo)) ® 5 a 0(-, a(oo))/ + S^i)/, 

^/iT(a(oo)) 

where the first operator is due to the resonance (and therefore must be rank one) 
and ||iSoo(i)/l|oo < (t) -1 ||/||wi,i. If we choose / to be a Schwartz function, say, 
then the first operator here is well-defined (despite the fact that d a <t> $ L 2 ) and it 
is given explicitly as 

c d a <f)(x, a) J d a (f>(y,a)f(y)dy. 

Let us give a naive reason for the constant term which the resonance produces in 
the sine-evolution. First, the sine-evolution operator acting on V0 (a zero eigen- 
function) would give a growth rate of t (just Taylor-expand). The dispersive part 
decays like t . Naively, the resonance should be half-way in between and grow 
something like t°. This is actually the case. In our proof, we show also that the 
cosine evolution operator is dispersive. 

We continue our discussion of the stable manifold theorem for NLW. As men- 
tioned previously, we make the ansatz 

(3) il>{x, t) = <p{x, a(t)) + u(x, t) 

The theorem is proved by solving for a and u globally in time in such a way that u 
disperses and a{t) converges to a limit which is 5-close to a(0) = 1. Note, however, 
that we cannot prove energy scattering for u to a free wave, which explains why we 
stated the scattering for a different decomposition of u. To understand this better, 
one checks from the Duhamel formula 1 for the nonlinear solution ip that for any 
t 2 > ti > 

^(•,i 2 )-VMi) £i 2 

Since 

^(.,ta) - i>{-M) = <f>(;h) - <j>(;h) + «(-,*a) - 
and (j>{-,t2) — (f>(-,ti) L 2 we see that 

u{-M)~u{-M)^L 2 

On the other hand, if u were to scatter to a free wave, then dtu G L 2 and conse- 
quently u(-, t-2) — u(-, ti) G L 2 . The truth is that u scatters to a free energy solution 
up to a small, decaying piece which has the same behavior at spatial infinity as 
d a (j)(-,a). In fact, our proof shows that 

(u,d t u) = (u,d t u) + (0,-a(t)d a <f>(;a(t))) +o^i xL2 (l) 

as t — > oo where (u,dtu) is a free wave with (u,dtu)(Q) 6 H 1 x L 2 . Consequently, 
if we set 

v(-,t) = u(-,t) + a(t)) - <f>(; o(oo)) 

then v scatters and 

i){-,t) = <t>{;a{oc)) +V(; t) 

as claimed in our theorem. The ansatz @ modifies the profile of <f> at spatial infinity 
instantaneously. While this is convenient technically, it may be less so physically, 
since perturbations should propagate at finite speed. The fact that we modify <fi 



^This was pointed out by Herbert Koch 
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outside of the light-cone of course forces us to absorb the change in that region 
by means of u. This is exactly why u does not scatter. It is natural to ask for 
an ansatz which does modify the profile of <fi only inside the light-cone. Thus, the 
question is as follows: can we write the stable solutions as 



where xi is a cut-off inside the light cone and \2 is a cut-off outside the light-cone, 
and u disperses and scatters in the energy space? It remains to be seen if our proof 
can accommodate this ansatz. 

One of our goals for this lecture is to explain how to obtain the ODE for a(t). 
This hinges on the linear estimates that we derived for the wave equation with a 
potential. During the past twelve years there has been much activity on dispersive 
and Strichartz estimates for the wave equation with a potential. Specifically, we 
would like to mention the work by Beals, Strauss, Cuccagna, Georgiev, d'Ancona, 
Pierfelice, Visciglia, and Yajima. However, these authors either assume that the 
potential is positive, or small, or that zero energy is neither an eigenvalue nor 
a resonance. In addition, varying degrees of smoothness and decay are assumed. 
Since we are dealing with a negative potential that leads to a zero energy resonance, 
none of these results apply to our problem. 

However, for the Schrodingcr equation with a potential there has been recent 
work of a similar flavor to what is needed here. Let us mention the following result 
from 2004 (there is a matrix version of this result from 2005 which is relevant to 
the linearized NLS, see above). 

THEOREM 5 

( Erdogan-Schlag ) 

For any Hamiltonian H = —A + V with a real-valued potential V(x) decaying like 
(x)^ 13 , (3 > 12, we have 



where F t is an operator for which sup t ||i 7 ' t ||x,i_ + L 00 < °°- 

Furthermore, if zero is not an eigenvalue but only a resonance, then Ft is rank 
one, and one can give an explicit expression for Ft. Independently, Yajima proved 
a similar result which required less decay of V and he gave an explicit expression 
for F t in all cases. 

As mentioned before, a unifying approach to resonances and/or eigenvalues at 
zero energy is given by the Laurent expansion of the perturbed resolvent around 
zero. To this end, we again recall the Jensen, Kato formula 



t) = <t>(;a(t)) X l(; t) + Cf>(; 1) X2 (-, t) + U(; t) 



||(e^P c -t-M)/||oo<H||/||i, 



(-A + V~z 2 ) 



2\-l 




-I h C + C\Z + ... 



z 



where we can incorporate the non-singular terms into 0(1) relative to a suitable 
weighted L 2 operator norm. To understand the relevance of this expansion, let us 
consider the free case V = in dimensions one and three. 
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Example 1 

Let d = 1 and V = 0. We then have, for all Imz > 0, 

]2 iz\x-y\ 

1 e "|a!-w| 



2«£ 2iz 

= ^ + 0(D. 

Hence, c_i = Jjl <g) 1. This tells us that 1 is a resonance function. Indeed, it is 
clear that y(x) = 1 is a solution of —y"(x) = 0. More generally, we say that zero is 
a resonance of + V provided there is a bounded, nonzero, solution of 

The reason for L°° lies with the mapping properties of the free resolvent. Also, if 
V decays sufficiently fast, then zero energy is never an eigenvalue (this is specific 
to d = 1). 

Example 2 

For d = 3 and V = 0, we have for Im z > 

p iz\x—y\ 

(-A- 2 2 )- 1 (^y) = fi r 

47r|a; — j/| 

Hence, we see that there is no singular term in z and hence no resonance at 0. As 
in the case d = 1 there is a characterization of resonances as solutions of Hf = 0. 
Indeed, we say that zero is a resonance iff there is / G L 2 '~ a \ L 2 for er > | 
(—A + V")/ = 0. In particular, /(a;) — 1 is not a resonance in <i = 3. To understand 
this definition better, let us write (—A+V)f = cquivalcntly as /+(— A) _1 V"/ = 0. 
Explicitly, 



m = -^[ P^-J{y)dy = ~f V(y)f(y) 
4?r J r3 \x-y\ 4tt J r3 



dy + 0((x)- 2 ). 



This shows us two things: 

• Since | a: | 1 : L 2 - ai — > L 2 ' _Q2 (R 3 ) as convolution operator provided ai, «2 > 
1/2 and ai + 02 > 1, we see that necessarily 

/eL v '(K 3 ) V,J> ^ 

at least if |I^(a;)| < (x)~ 2 ~ e . This explains the definition of a resonance 
function. 

• In fact, in that case / decays like at infinity iff J Vf(y) dy 7^ (the 
resonant case), whereas / is an eigenfunction decaying like \x\~ 2 at infinity 
iSjVf(y)dy = 0. 

It is of course important to understand how to obtain Laurent expansions of the 
perturbed resolvent. This was done by Jensen and Kato around 1978, but a more 
transparent approach was recently found by Jensen and Ncnciu. Let us state the 
main lemma from their paper. 
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LEMMA 3 

( Jensen-Nenciu - 2001 ) 

Let H be a Hilbert space, and suppose that bounded operators 

A(z) = A + zA 1 (z) :H^H 

are given for all z £ F C C \ {0}. In addition, assume that A\(z) is uniformly 
bounded in z, and that Aq = Aq. Suppose that zero is an isolated point of the 
spectrum of Aq and let 

S = Proj(KerAo) 

be of finite, positive rank. Then (A Q + S*) -1 exists, and thus (A(z) + S)^ 1 also exists 
for small z. Define 

B(z) = -(S-S(A(z) + S)- 1 S) 
z 

Then A(z) is invertible iff B(z) is invertible on SL 2 , and in that case we have 
A(z)- 1 = (A(z) + S)- 1 + -(A(z) + S , )- 1 55(z)- 1 S(A(z) + S)' 1 . 

Note that B(z) is uniformly bounded in the operator norm as z — > 0. Indeed, 

S - S(A Q + s)-^ = 

due to the self-adjointness of Aq. 

The connection between this lemma and resolvent expansions is furnished by the 
symmetric resolvent identity: Let V = v 2 U were U = sgnV. Then 

R v {z) := (-A + V - z 2 )- 1 = R (z) - R (z)v(U + «i? (z) V )- 1 U i? (^) 

for all Imz > 0. Let us set 

A(z) = U + vR, (z)v = U + vRo(0)v + z <Mz) ' Ro(0))v = _ ^ + 

z 

If Aq is invertible, then the Laurent expansion of Ry{z) around z = has no 
singular terms. The converse is also true. Thus zero energy is neither an eigenvalue 
nor a resonance iff Aq is invertible. Using Fredholm theory (and the compactness of 
vRq(Q)v on L 2 ), we conclude that Aq is invertible iff there is no nonzero L 2 solution 
of 

(U + v(-A)- 1 v)f = 

This in turn is equivalent to the existence of certain solutions h to 

(-A + V)h = 

Indeed, formally set h = (— A)~ 1 vf. 

In our application to NLW, Aq is therefore not invertible. We apply the Jensen, 
Nenciu lemma and check that the B{z) is invertible. Since S is rank-one in this 
case, this invertibility is equivalent with the nonvanishing of a number. Not sur- 
prisingly, this number exactly turns out to be J V(y)d a (f>(y, a) dy ^ 0, see above. 
Hence, the Jensen, Nenciu lemma allows us to perform the Laurent expansion with 
relatively explicit coefficients. This latter feature is most relevant for understanding 
the evolution operators of the perturbed wave equation. Indeed, by the spectral 
calculus, 

(4) ^f D = ir ^R v (X)XdX. 
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In order to derive dispersive estimates we follow the usual path of splitting the 
integration here into |A| > e and |A| < e for some e small enough. It is common 
knowledge that dispersive estimates for perturbed operators belong to the realm 
of a zero energy theory. The point here is that energies different from zero cannot 
destroy the free rate of decay. Only zero energy can do that. Nevertheless, the 
regime |A| > e is non-trivial in particular with respect to the amount of decay of V 
required to make this work. Recall that we have exactly \x\~ 4 decay in NLW. To 
appreciate the difficulty of lowering the decay requirements on V see Goldberg's 
paper |Gol| . At any rate, in this regime we work with the (finite) Born expansion 



Plugging the terms here not involving Ry into the right-hand side of (J3J leads 
to completely explicit expressions for the kernels due to the explicit nature of Ro 
in R 3 . It is therefore possible, with some work, to obtain the desired dispersive 
bounds on those kernels. For the remainder term, which still involves Ry, we rely 
on the so-called limiting absorption principle due to Agmon. This refers to the fact 
that the (perturbed) resolvent remains bounded on weighted L 2 spaces. For further 
details on not small energies, we refer the reader to our paper. 

For energies |A| < e, we expand the resolvent Ry as in the Jensen, Nenciu lemma. 
The singular term ultimately leads to the non-decaying part 



whereas we show in our paper that the nonsingular terms yield a dispersive operator 
S(t), see (J2J above. Dispersion here refers to bounds in L°°. It is an open question 
whether or not S(t) satisfies Strichartz estimates. For this see the author's work 
on Littlcwood-Palcy theory and the distorted Fourier transform. 

To conclude these lectures, let us turn to the question of how to determine the 
ODE for a(t). To this end let us recall that we write the solution vector in the form 



As explained before, generally speaking n+ is an exponentially unstable mode. It 
needs to be stabilized by means of the stability condition This is the origin of 
the codimension one condition. On the other hand, n_ is exponentially decaying. 
Finally, the remainder u satisfies the equation 

dttu + H(a(oo))u = 



(5) [-d tt <t>(; a(t)) + (V(-,a(oc)) - V(; a(t))u(- : t) + N(u(- : t),<j>{; a(t))} 



It is our goal to show that u (and therefore also u) disperses, but this is non-obvious 
because of the non-dispersive term in (J2J). This is precisely where a(t) comes in, 
the point being that we will collect all non-decaying terms and observe that they 
are all multiples of ip := d a (f>{- , a(oo)) . To give a flavor of this, let us consider 
the contribution of the first term in J^J to the Duhamel formula. It is (writing H 
instead of H(a(oo)) for simplicity) 



Rv =Rq- RqVRo + RoVRqVRq t •■■ ± RaV • • • VR V V ■ ■ ■ VR 



cod a cj)(-,a) (g) d a <i>{-,a) 




I 



* sin((t - s)VH) 

7h 



P g<x> d s a{s)d a (j>(-,a(s))ds 
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We could use the representation J5J here. However, this is very misguided since it 
is impossible to let 

C 1p (g) "0 

act on d a 4>{- , a(s)) (because of d a <fi $ L 2 (M. 3 )). Instead, we proceed as follows: 







ri V ti 



(6) + / d(s)cos((t-s)Vff)P^A^(-> a ( s )) ds 

Now observe that 



i] 



Using these relations and the fact that (for a±, a 2 € (1/2, 2)) 

d a <f>(;ai) = {a 2 /ai)i d a ^{-,a 2 ) + 0(K - a 2 |(a;)" 3 ) 

we reduce d a (j>(-, a(s)) to tp by means of the scalar factor (a(oo)/a(s))*. This shows 
that we need d(0) = 0, and also allows us to pull out ip from the integral above with 
a factor depending on a(t). The remainder which comes from the 0((x) _3 )-term 
turns out to be dispersive. We refer the reader to our paper for further details. 
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